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Abstract 

The radiative ip — > tt^tjj decay is discussed emphasizing the effects 
of the ao(980) scalar resonance which dominates the high values of 
' the TT^r] invariant mass spectrum. In its lowest part, the proposed 

amplitude coincides with the reliable and ChPT-inspired contribution 
coming from chiral loops. The ao(980) resonance is then incorporated 
exploiting the complementarity between ChPT and the linear sigma 
model for this channel. The recently reported experimental invariant 
mass distribution and branching ratio can be satisfactorily accommo- 
dated in our framework. For the latter, a value of B{(j) tt^t]^) in 
the range (0.75-0.95) X 10"^ is predicted. 



1 Introduction 



The 1 GeV energy region is a particularly challenging domain. On one side, 
it is far below from the perturbative QCD regime and, on the other, Chiral 
Perturbation Theory (ChPT) is not expected to make reliable predictions at 
these energy values where resonance effects are known to be present. Among 
the latter, those proceeding by the exchange of the scalar resonances /o(980) 
and ao(980) should dominate their respective channels. The controversial 
nature of these resonances and the poor knowledge of their properties 
[0] adds then further complexity (and interest) to this 1 GeV energy region. 
Indeed, several proposals have been suggested along the years concerning the 
constitution of these scalars as complex qqqq states ||^, KK molecules [0 or 
ordinary qq mesons f^. 

The Novosibirsk CMD-2 and SND Collaborations have reported very re- 
cently, among others, the branching ratio and the tt^t] invariant mass dis- 
tribution for the — > TT^rj'-f decay. For the branching ratio, the CMD-2 
Collaboration reports 5(0 ^ TT°r]-f) = (0.90 ± 0.24 ± 0.10) x 10"^ [§, while 
the SND result is, consistently, B{(f) tt^t]-^) = (0.88±0.14±0.09) x 10"^ [§. 
The observed invariant mass distribution shows a significant enhancement at 
large n^r] invariant mass that, according to Refs. ||^, 0], could be interpreted 
as a manifestation of a sizeable contribution of the ao (980)7 intermediate 
state. This and other radiative decays are also expected to be intensively 
investigated at the Frascati 0-factory DA$NE 

On the theoretical side, the V p^p^'y decays have been considered 
by a number of authors |^, |TU|, [TT|, |T^. In particular, it has been shown 
that the intermediate vector meson contributions to — ir^rj'j lead to a small 
_B(0 — i> 7r°r77)vMD = 5.4 x 10^^ [|13|, whereas a chiral loop model closely linked 
to standard ChPT predicts 5(0 'r^^Vl)x = 3-0 x 10"^ |14l. Needless to say. 



the scalar resonance effects and, in particular, the resonance pole associated 
to the ao(980) were not contemplated in these two approaches. The recent 
experimental data from Novosibirsk — for both the branching ratio and the 
TT^Tj invariant mass spectrum showing an enhancement around the ao(980) 
mass — seem therefore to disfavour these predictions based on vector meson 
exchange and/or a simple extrapolation of ChPT ideas. 

If we rely on the resonance picture, it is clear that the ao(980) scalar 
meson — lying just below the mass and having the appropriate quantum 
numbers — should play an important role in the 7r°?77 decay. Several 
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theoretical attempts to describe the effects of scalars in radiative decays 
have appeared so far. Among others, we would like to refer to the "no 
structure" model |T3, to the K^K" model pi, 0, where the — > ao7 



amplitude is generated through a loop of charged kaons, and to the chiral 
unitary approach (UxPT) where the decay — > ir^'q'y occurs through 
a loop of charged kaons that subsequently annihilate into TT^rij. In the two 
former cases the scalar resonances are included ad hoc while in the latter 
they are generated dynamically by unitarizing the one-loop amplitudes. 

In this letter, we are mainly interested in incorporating scalar resonances 
and their pole effects into a ChPT inspired context [|16]. While vector 
and axial-vector resonances can be included in a transparent and success- 
ful way, offering some theoretical basis to conventional vector meson dom- 



inance (VMD) ideas ||T^, the incorporation of scalar resonances has been 
more ambiguous and less successful up to now []TB| . In order to take explic- 
itly into account scalar resonances and their pole effects, we propose to use 
the linear sigma model (LcrM). This will allow us to take advantage of the 
common origin of ChPT and the LcxM to improve the chiral loop predictions 
for V P^P^'j exploiting the complementarity of both approaches for these 
specific processes. On one side, ChPT is the established theory of the pseu- 
doscalar interactions at low energy. However, it is not reliable at energies of 
a typical vector meson mass and, as just stated, scalar resonance poles are 
not explicitly included. As a consequence, ChPT inspired loop models can 
give rough estimates for B{V P^P^'j) but will hardly be able to repro- 
duce the observed enhancements in the invariant mass spectra. On the other 
side, the LcrM is a much simpler model dealing similarly with pseudoscalar 
interactions but incorporating scalar resonances in a systematic and definite 
way. Thanks to this, the LaM should be able to reproduce the resonance 
peaks in the spectra and, although it does not provide a systematical frame- 
work for the pseudoscalar meson physics, this model could be of relevance 
in describing the scalar resonances when linked to a well established ChPT 
context. In order to show in detail the proposed framework, we will focus our 
attention on the — >■ 7r°?77 decay mode. Other decay modes are somewhat 
more involved and will be analyzed in forthcoming work. 
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2 (p ^ TT^r]^ and chiral loops 

The vector meson initiated V p^p^'y decays cannot be treated in strict 
Chiral Perturbation Theory (ChPT). This theory has to be extended to in- 
corporate on-shell vector meson fields. At lowest order, this may be easily 
achieved by means of the ChPT Lagrangian: 

C2 = ^{D.U^D'^U + U^x + X^U) , (1) 

where f = = 92.4 MeV at this order, U = exp{i\/2P/ f ) with P being the 
usual pseudoscalar nonet matrix, and x = 2BqM. with A4 = diag(mu, nid, rris). 
The covariant derivative, now enlarged to include vector mesons, is defined 
as D^U = d^U - ieA^[Q,U] - ig[V^,Ul with Q = diag(2/3, -1/3, -1/3) 
being the quark charge matrix and the additional matrix containing the 
nonet of ideally mixed vector meson fields. The diagonal elements of V are 
(p° + uj)/\^, (— p° + uj)/\/2 and 0, thus following the same conventional 
normalization as for the pseudoscalar nonet matrix P. 

There is no tree-level contribution from this Lagrangian to the — > ir^rj'y 
amplitude and at the one-loop level one needs to compute the set of diagrams 
shown in Ref. |]14|. A straightforward calculation leads to the following finite 
amplitude for 0(g*,e*) 7r°(p)r7(p')7(g, e) (see Ref. [0 for further details): 

-A(0^A7)x = ^^2fr^(e*eg*g-e*geg*)i:(m2o^) 

^+ (2) 



X A{K+K' ir^ri) 



where e*t q*q — t*q tq* makes the amplitude Lorentz- and gauge-invariant, 
m^o^ = s = {p+p'Y = {1* is the invariant mass of the final pseudoscalar 
system and L{m\o ) is the loop integral function defined as 
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The coupling constant 
) = ge*{p^ — pJ) with 



g comes from the strong amphtude A^cj) 
\g\ = 4.59 to agree with 1(0 K+K~)^^p = 2.19 MeV. The latter is the 
part beyond standard ChPT which we have fixed phenomenologically. The 
four-pseudoscalar amplitude is instead a standard ChPT amplitude^ which 
is found to depend linearly on the variable s 
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(5) 



In the calculation of the decay amplitudes (Q) and (||) we have introduced 
ri-Vj mixing effects. As it is well known, a rigorous and general extension of 
5*^/(3) ChPT to include the ninth pseudoscalar meson r/o is not straightfor- 
ward and requires the introduction of new terms in the chiral Lagrangian 
|l^. However, if one relies on classical arguments based on nonet sym- 
metry, a phenomenologically successful description of the r/-//' system is 
achieved [|1^]. The 77-17' mixing angle is then found to be compatible with 
Qp = arcsin(— 1/3) ~ —19.5°, quite in agreement with recent phenomenolog- 
ical estimates |2^. In Sect. ^, it will be shown that this choice for the 77-77' 
mixing angle greatly simplifies the calculation of the — > Ti^rj'y amplitude in 
the LctM, and reduces up to a minimum the number of free-parameters. 
The invariant mass distribution for the —* 7r°7/7 decay is predicted to 



1 A{K+K- ^ TT^rjs)^ 
account as in Ref. [[l4|. 
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Figure 1: dB{(f) 7r°r77)/(im^o^ x 10^ MeV ^ as a function of the mT^o^ 
invariant mass in a chiral loop model. Experimental data are taken from 
Ref. m. 



be given by the following spectrum (see Fig. |1|): 

(6) 

Integrating Eq. (P) over the whole physical region one obtains for the branch- 
ing ratio: 

B{(P A7)x = 0-47 X 10-^ . (7) 

As expected, Fig. ^ shows that our chiral loop approach gives a reason- 
able prediction for the lower part of the spectrum but fails to reproduce the 
observed enhancement in its higher part, where ao(980)-resonance effects (ig- 
nored up to this point of our approach) should manifest. As a consequence, 
the predicted branching ratio turns out to be below the experimental value 
by about a factor of 2. 
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3 Improved approach to ^ n^r]^ 

To analyze the scalar resonance effects in the V — > P^P^'-f decay amplitudes, 
the linear sigma model (LcrM) |2l| will be shown to be particularly appro- 
priate. It is a well-defined U{3) x U{3) chiral model which incorporates ab 
initio both the nonet of pseudoscalar mesons together with its chiral part- 
ner, the scalar mesons nonet. Recently, the model has been resurrected as a 
framework to study the implications of chiral symmetry for the controversial 
scalar sector of QCD, and some variations of the basic LaM Lagrangian have 
shown to be phenomenologically rather successful |2^, |23|, 



In this context, the V — > P^P^j decays proceed through a loop of charged 
pseudoscalar mesons emitted by the initial vector. Due to the additional 
emission of a photon, these charged pseudoscalar pairs with the initial J^'-' = 
1 quantum numbers can rescatter into J^'" = pairs of charged or 
neutral pseudoscalars. The scalar resonances are expected to play an essential 
role in this rescattering process and the LcxM seems mostly appropriate to 
fix the corresponding amplitudes. 

Several simplifications happen when one considers the TT^rj'y decay 
mode. As in the analysis of Sect. |^, contributions from charged pions in the 
loops are highly suppressed because they involve the isospin violating and 
OZI-rule forbidden ^vrvr coupling; hence, the dominant contributions arise 
exclusively from loops of charged kaons. The subsequent rescattering of these 
charged kaon pairs into the final J^*" = 0^^ n^r] state is then quite simple. 
Indeed, the LcrM amplitude for K^K^ -k^t] contains a contact term, a 
term with an oq exchanged in the s-channel and two terms with a k {i.e. the 
strange 7 = 1/2 scalar resonance) exchanged in the t- and w-channels. How- 
ever, the latter /t-exchange contributions are absent for an rj-rj' mixing angle 
Op = arcsin(— 1/3) ~ —19.5° since the g^Rrj coupling constant appearing in 
one of the k vertices vanishes. The calculation is then reduced to the dia- 
grams shown in Fig. ^, and it is thus much simpler for this particular (but 
phenomenologically acceptable ||20|) choice of the rj-rj' mixing angle. Moder- 
ate departures from this value translate into weak Qt^Kri couplings^ appearing 
in the — > TT^rj'-f amplitude. Their effects seem to be small (see below) for 
the present and expected levels of experimental accuracy and for this kind 



^The coupling constant g^Kr] is proportional to sin [0p — arcsin(— 1/3)] and its depen- 
dence on the rj-rj' mixing angle for values around 9p = arcsin(— 1/3) ~ —19.5° is soft. 
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Figure 2: One- loop Feynman diagrams for (j) n^rj'j in the LaM for a rj-rj' 
mixing angle dp = arcsin(— 1/3) ~ —19.5°. 

of processes governed by poorly known scalar resonances. More importantly, 
the absence of k contributions makes our predictions more definite and solid 
since we avoid one of the major uncertainties affecting the scalar nonet dy- 
namics, namely, the mass of its strange members. Indeed, recent analysis by 



various authors require a light fi:(900) [0, |2^, ^ while other authors 
deny the existence of a low mass pole |^ and identify the PDG ii'Q(1430) 
state with the strange member of the lowest lying scalar nonet ||2^, |2^ . 

A straightforward calculation of the (p ^ ir^rj'-f decay amplitude leads 
to an expression identical to that in Eq. but with the four-pseudoscalar 
amplitude now computed in a LaM context. In this case, the amplitude is 
just 

A{K^K- vr%)L.M = gK^K-.o, - ^""^"f^f" , (8) 
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where -Da„(m^o^) is the propagator and the couphng constants areQ 

iJK+K TT^rj 9f,^ Y 3 2fKf-n ' 



9aoK+K 2/k ■ 

This amphtude can then be rewritten as 

A{K^K~ - A)l.m = ^^(^^0, - ml) X . (10) 



We would hke to make a few remarks on the four-pseudoscalar amphtude 
in Eq. (|10D and compare it with the ChPT amphtude in Eq. (^): 

i) For rria,, — > oo and ignoring SU (3)-breaking in the pseudoscalar masses 
and decay constants, the LaM amphtude ( [T(]| ) reduces to the ChPT 
one (^. The former consists of a constant four-pseudoscalar vertex 
plus a second term whose s-dependence is generated by the s = m?^o^ 
piece in the ao propagator Da^is = ""^^o^), as shown in Eq. (||). Their 
sum (see Eq. ([To|) ) in the good SU{3) and niao ^ limits ends up 
with an amplitude which is linear in s = m'^o^ and mimics perfectly the 
effects of the derivative and massive terms in the ChPT Lagrangian (|^) 
leading respectively to the two terms in the ChPT amplitude (^). This, 
we believe, is the main virtue of our approach and makes the use of the 
LaM reliable at least for amplitudes like ours where s-channel exchange 
plays the main role and {t, M)-channel exchange can be ignored. 

ii) The LcxM and ChPT yield slightly different amplitudes in the m^y oo 
limit because of the way 5'?7(3)-symmetry is broken in the two ap- 



proaches. In the case of the LaM |22, 23, 24|, a non SU{3) sym- 
metric choice of the vacuum expectation values makes simultaneously 
7^ and fn 7^ fx, whereas in ChPT 7^ is already present 
in the lowest order Lagrangian while /tt 7^ fx is only achieved at higher 
orders. 

iii) In the (p Tc^rij decay, the threshold for tt^t] production is not far 
from the mass of the ao(980) and that makes crucial the incorporation 



^See Ref. [p2| for a detailed calculation of these coupling constants in the LaM. 
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of the ao in an explicit way. Due to the presence of the full propagator 
DaQ^s), as in Eq. ([lOD , such an amplitude — closely linked to that from 
ChPT and thus expected to be able to account for the lowest part of 
the n'^T] mass spectrum — should also be able to reproduce the effects 
of the ao pole at higher n^r] invariant mass values. 

iv) The need for the propagator introduces, however, some uncertain- 
ties in our treatment. Indeed, the opening of the KK channel near the 
ao(980) mass has motivated the use of different expressions for /^^^(s). 
A first possibility consists in using a Breit-Wigner propagator with an 
energy dependent width (to incorporate the known kinematic correc- 
tions) : 



where 



m 



ao 



(11) 
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9{^/s - 2mK+] 



+ 



1 - 
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6'(yi — 2m^o) 



(12) 

Another interesting and widely accepted option was proposed by Flatte 
time ago specifically for the two-channel resonance pOf. The relative 
narrowness of the observed vri] peak around 980 MeV is then explained 
by the action of unitarity and analyticity at the KK threshold. This 
amounts to extend the preceding formulae below the KK threshold to 
include the now purely imaginary kaon contributions. 

Due to these distinct possibilities to deal with the propagator, as well 
as to other differences introduced when implementing and fitting the basic 
LcrM Lagrangian by several authors, a set of predictions can be obtained 
for the four-pseudoscalar amplitude (|10D . In turn, these various amplitudes 
have to substitute the four-pseudoscalar ChPT amplitude in Eq. (D) to finally 
obtain the corresponding invariant mass distributions of the (j) — > 71*^777 decay 
modeFl Our purpose hereafter is to briefly discuss a few of these treatments 



*Since all the four-pseudoscalar amplitudes we are considering depend only on the 
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in order to show that the observed properties for this specific decay can be 
accommodated in our ChPT- and LcrM-inspired approach. 

We start our discussion along the hues of Ref. taking for the Oq 
propagator the simple Breit-Wigner prescription in Eq. (|lT]). The use of 



this propagator for the LcrM amplitude Eq. ([T0|) and its insertion in Eq. 
predicts the m^o^ invariant mass spectrum shown by the dotted line in Fig. ^ 
Integrating over the whole physical region leads to the branching ratio 

B{<P ^ 7rV)L<xM[22] = 0.80 X 10-^ , (13) 

in agreement with the experimental data ||^, |^. However, since the simple 
expression used for the propagator implies a large oo-width (Fa^^^^ ~ 460 
MeV [^), the desired enhancement in the invariant mass spectrum appears 
in its central part rather than around the Oq peak. 

This unpleasant feature is easily corrected when turning to the proposal 
by Tornqvist Indeed, a Gaussian form factor related to the finite size 
of physical mesons and depending on the final CM-momentum is introduced 
to describe the decays of scalar resonances in this approach. As a result, the 
decay width of ao(980) into n^rj is reduced {Vao^-wr) — 273 MeV p3|) without 
affecting that of ao(980) into KK. This fact produces an enhancement in 
the spectrum for the higher values of the ^,^0^ invariant mass, as shown by 
the dashed line in Fig. ^. The integrated branching ratio is then predicted 
to be 

5(0 ^ vrV)L.M[23] = 0.90 X 10-^ , (14) 

in good agreement with the experimental data [§,0]. A possible difficulty 
of this approach at the phenomenological level is that it predicts an rj-r]' 
mixing angle of 9p ~ —5.0° [^], considerably less negative than the usually 
accepted value {6p between —20° and —15° and the value required in our 
simplified analysis 9p ~ —19.5°. This allows for an estimate of the typical 
errors introduced when neglecting K-exchange in the t- and -u-channels as 
compared to ao-exchange in s-channel. The significant factor is the ratio of 
coupling constants r = {gf,±K±rj gK±K±no) / {gaonOr, gaoK+K-), which vanishes 
in the ideal situation where 6p ^ —19.5° but takes the value r ~ — 1/3 for 
6p ~ —5.0° if one ignores S'?7(3)-breaking corrections. To the smallness of 
r ~ —1/3 one has to add the fact that the amplitude for the almost on-shell 

variable s, they factorize out of the loop integration and the structure of Eq. is fully 
preserved. 
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Figure 3: dB((f) — » n'^'q'j) j ^771,^0^ x 10'' MeV~^ as a function of the 7717^0^ invari- 
ant mass in the LcrM. The dotted, dashed and sohd hues correspond to the 
versions of the LcrM proposed by Refs. p2|, ^ respectively. Experimental 
data are taken from Ref. [[7|. 



ao-exchange in the s-channel is mainly imaginary and does not interfere with 
the almost real amplitude for off-shell ^-exchange in the t- and w-channels. As 
a result, the error in Eq. (|14D introduced by neglecting this latter contribution 
can be estimated to be below some 10% even for such unusual values of 
Qp - -5°. 

None of these drawbacks are encountered when turning to the treatment 



proposed by Shabahn ||2^. In the fitting procedure adopted by this author 
no attempt is made to fix the 77-77' mixing angle within the model. Thanks to 
this, one minimizes the uncertainties associated with the incorporation of the 
ninth pseudoscalar meson 770 via the axial anomaly term. The value of the 
77-7/' mixing angle is then fixed outside the model to its phenomenologically 
preferred value Qp ~ —19.5°. Another relevant feature of Shabalin's approach 
is the introduction of the well-known Platte corrections to the ao(980) prop- 
agator. The ao-width is then drastically reduced from the uncorrected value 
Y ao — 304 MeV to a more acceptable visible width of — 65 MeV. With 
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all this information taken from Ref. 
variant mass spectrum shown by the solid line in Fig. 
branching ratio 

5(0 ^ vr V)l<xM[241 = 0.93 x 10"^ 



our approach predicts the m^o^ in- 
and the integrated 

(15) 



Both the spectrum and the branching ratio are in nice agreement with the 
experimental data P, |^. The fact that Shabalin's model incorporates the 
Flatte corrections to the oq resonant shape has played a substantial role 
in this achievement. 



4 Conclusions 

The main aim of the present letter has been to propose and discuss an am- 
plitude for the radiative — » 7r°?77 decay exploiting the complementarity 
between ChPT and LaM ideas. Thanks to the latter, our amplitude con- 
tains the full propagator of the ao(980) scalar resonance which dominates 
the higher part of the ir^r] invariant mass spectrum. In the low invariant 
mass region, where ChPT is expected to work quite reliably, the proposed 
amplitude is shown to coincide with that coming from a chiral loop calcu- 
lation. This, we believe, makes reliable our approach to the V p^p^'j 
dynamics and, in particular, to the — > n^ri'-f decay mode for which some 
simplifying conditions hold and lead to a simple and well-defined amplitude. 
Then our predictions depend only on a reduced number of parameters which, 
in principle, can be extracted from independent data. Some of these data 
refer to scalar meson properties which are not well established and thus affect 
the accuracy of our predictions although by no means in a drastic way. 

We can safely conclude that all the reported properties for the — *■ 
7r°?77 decay mode can be accommodated in our approach. The branching 
ratio is predicted to be in the range BR{(f) n'^ri'y) = (0.75-0.95) x lO"'^, 
compatible with the present available data. Similarly, the measured ir^r) 
invariant mass spectrum is reproduced by our amplitude in a reasonable 
way. The uncertainties affecting these predictions suggest that further tests 
and more refined analyses are needed, particularly when the higher accuracy 
data from ongoing experiments will be available. This should contribute to 
clarify one of the most controversial aspects of hadron physics: the scalar 
states around 1 GeV. 
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